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1. Introduction 

The large- N expansion is an alternative to conventional perturbation theory. In both 
QCD and supersymmetric QCD (SQCD), the large N c expansion is of particular interest 
[l]; more tractable, however is the large Nf expansion, and it is this that we study here, 
for SQCD. In any theory, calculation correct to some non-trivial order in 1/N requires the 
summation of one or more infinite sets of Feynman diagrams, and hence possible insight 
into the large-order behaviour of perturbation theory. There is a considerable literature 
devoted to the study of large order perturbation theory 0, but remarkably little of it 
concerned with supersymmetric theories. [For a discussion of the large order behaviour of 
the supersymmetric anharmonic oscillator, see Ref. 0; and for investigation of the role 
of renormalons in various supersymmetric theories see Ref. |4]] . Supersymmetric a-models 
have been studied at large iV using critical methods f| in a series of papers by Gracey: 
see for example Ref. ||.] In this paper we calculate the leading and 0(1/Nf) terms in 
the gauge beta-function, g , for SQCD0. The type of bubble-sums we confront are in fact 
similar to those in Ref. Jl ; the difference being that they were concerned with renormalon 
singularities in amplitudes, whereas, our interest being in beta-functions, we require the 
ultra-violet divergent terms from such sums. We find that all our results can be expressed 
in terms of a simple function of the coupling constant, and that although infinite classes 
of diagrams have been summed, the resulting coefficient of 1/Nf has a finite radius of 
convergence in g. Of course we cannot expect the expansion in powers of 1/Nf to be 
convergent; it might well, however, be Borel-summable, as has been conjectured || for the 
expansion in g. 

Some of our results have already appeared in Ref. || ; here we give more calculational 
detail and also extract the SQCD case. We use the superfield formalism, allied with su- 
persymmetric dimensional regularisation and minimal subtraction (DRED)[]1(J. As always 
with superfield perturbation theory, the calculation of a typical contribution consists of first 
reducing a superfield Feynman diagram to a normal one by performing D-algebra, and then 
performing the resulting Feynman integral. The first part is straightforward, because our 
diagrams consist of bubble chains grafted on to otherwise (at most) three-loop graphs. The 
second part appears formidable; but fortunately the key to its performance has been pro- 
vided by Palanques-Mestre and Pascual[[ll|, who carried out similar non-super symmetric 
calculations in the abelian case. The crucial realisation is that the automatic cancellation 

1 For QCD at large Nf see § 
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of non-local counter-terms leads to apparently miraculous (but easily verified) identities 
which simplify the summation over subtractions. 

In the next section we carry out our program for an Abelian theory, with a super- 
potential selected so that Af = 2 supersymmetry is included as a special case. Then in 
section 3 we show how with very little extra work we can extract the corresponding result 
for the non- Abelian case, by exploiting the fact that M = 2 theories are finite beyond one 
loop. In section 4 we show how the SQCD result can be also deduced, and in section 5 we 
compare our result with perturbation theory in g for the SQCD case. 



2. General Abelian Theory 

Here we present results for a general theory produced by a U\ gauging of the Wess- 
Zumino model with a superpotential given by 

A Nf 

W=^=J2^iXi- (2-1) 



v N f • , 

V J i = l 

Suppose we require a gauging such that the charges q Xi = q x and q^ = q^ are independent 
of i. Charge conservation and anomaly cancellation then require 

+ + = Ql + Nf(q 3 x + ql) = 0, (2.2) 

yielding 

(N f - l)(gj + ql) - 3(<gg c + q x ^) = 0. (2.3) 

An obvious solution is q% = —q x , q^ = 0. In fact if this does not hold, then dividing 
Eq. (|2.3|) by q x + q^ we obtain 



(JV> - l)(q 2 x + ql) - (N f - 2) M = 0, (2.4) 

which is easily seen to have no solutions (except in the case Nf = 1, with either q^ = 
or q x — 0, which are clearly not interesting from our present large- Nf perspective). So we 
now set q^ = and q% = —q x = q; we will take q = gj y/Nf. For the special case A = \f2g 
we have M = 2 supersymmetry. In Fig. 1 we show the Feynman diagrams we require for 

Ar 
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Fig. 1: Feynman diagrams for section 2. Wavy lines are vector propa- 
gators, solid lines are £ or x propagators, and dashed lines are <p prop- 
agators. Blobs denote sums of chains of£,x bubbles. 

A blob on a diagram represents the sum of chains of £, x bubbles of arbitrary length, 
as shown in Fig. 2(a) (for a (^-propagator) and Fig. 2(b) (for a vector propagator). We 
choose to work in the Landau gauge; this means that the diagram with no bubble on a 
vector propagator can be subsumed with the diagrams with chains of bubbles, since a chain 
of one or more bubbles produces a tranverse projection operator. Fig. 1(a) is 0(1), while 
Figs. l(b)-(g) are all 0(1/JV». 

-Q--0 O O 

(a) 




(b) 



Fig. 2: Bubble chains on chiral and vector propagators. 

For details of our technique for dealing with the D-algebra part of the calculation 
we refer the reader to Ref. []12], where we found (3 g for an abelian theory to four loops 
by calculating the vector superfield self-energy. The upshot is that each diagram with 
a blob reduces to one of a basic set of Feynman integrals, depicted in Fig. 3. Again in 
Fig. 3, a chain of bubbles is understood to represent a summation over arbitrary numbers 
of bubbles. 



4 




Fig. 3: Feynman diagrams representing the bubble sums A,B C and D. 
The black dots denote squared propagators and pairs of arrows denote 
contracted momenta. 

The Feynman integrals (A)-(D) are evaluated in Appendix A, with results given in 
Eqs. (A. 17), (A. 8), (A. 18) and (A. 25) respectively. We note all the bubble sums relevant to 
our calculations depend on the function G(x), which has a zero at x = 1 and a simple pole 
at x = |. We may therefore anticipate that our results in this and subsequent sections will 
have a finite radius of convergence in the appropriate coupling constant, because of this 
pole. The contributions from the diagrams in Fig. 1 are given in terms of the diagrams of 
Fig. 3 in Table 1. In this Table and the following ones, the second column represents the 
result of performing the D-algebra, and should be multiplied by the corresponding factor in 
the third column which represents the symmetry factor and also (in the non-abelian case) 
products of group matrices. An additional factor of l/Nf is also understood in each case, 
and we have introduced K — g 2 /(8TT 2 ) and y = A 2 /(167r 2 ). For Af = 2 supersymmetry we 
have y = K. 



Diagram 


Bubble Sum 


Factor 


1(b) 


\C{K) - 2A(K) 


2gK 2 


1(c) 


\[±B{K) - A(K)} 


4gK 2 


1(d) 


\[±B{K) - A{K)\ 


-8gK 2 


1(e) 




2gK 2 


1(f) 


A(y) 


2gKy 


1(g) 


C{y) 


-gKy 



Table 1: 0(1/Nf) contributions to j3 g for the abelian theory. 



Figs. l(b)-l(g) evidently have precise counterparts in a three-loop calculation with 
the blob representing just one bubble; and in fact the bubble-summed diagrams (A)-(D) of 
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Fig. 3 are generalised forms of those in Fig 2(A)-(D) of Ref. [12]. The results in Table 1 for 
Figs. l(b)-l(g) can hence be read off from the corresponding results of Table 1 in Ref. |T3 . 
Using our results from Appendix A for Fig. 3(A)-(D) and including the 0(1) contribution 
of Fig. 1(a), we obtain 



P g = gK 



2 r K 

1 + — J (1 - 2x)G(x) dx 



The anomalous dimensions through 0(1/Nf) are given by Fig. 4. 




(2.5) 



(a) 



(b) 



(c) 






Id) 



(e) 



(f) 




Ig) 



Fig. 4 : Feynman diagrams for the anomalous dimensions 

The contributions from these diagrams are given in Table 2, in terms of the basic 
diagrams of Fig. 3. 



Diagram 


Bubble Sum 


Factor 


4(a) 


B(y) 


y 


4(b) 


B{K) 


-K 


4(c) 







4(e) 


My) 


2y 2 


4(f) 


±[B(K)-1]-A[K) 


2yK 


4(g) 




-2yK 



Table 2: 0(1/ Nf) contributions to anomalous dimensions for the abelian 
theory 

The contributions from graphs with one or more bubbles in Fig. 4(c) are guaranteed to 
be zero since the bubbles lead to a transverse projection operator on the vector propagator, 
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giving zero when attached to an external line. Since we have chosen the Landau gauge, 
the zero-bubble contribution to Fig. 4(c) will also give zero by the same token. The 
same reasoning will be used later to deduce that Figs. 6(b), (d) give zero. Adding the 
contributions from Figs. 4(a)-(c), and those from Figs. 4(d)-(g), we obtain 



= 7x 



7^ 



i [yG(y) - KG(K)) 



y+ w f 



G(K) - G(y) + 2 f G(x) dx 

J V 



(2.6) 



It is easy to verify that our result for (3 g in Eq. (|2.5|) reproduces the relevant terms 
in the three and four loop calculations presented in Ref. []T2|. The results for x and 7^ 
in Eq. ( pip agree with the three- loop results of Ref. [TB| , and in the ungauged case agree 
with the four-loop results of Ref. for a generalised Wess-Zumino model. Moreover, 
for Af = 2 (which corresponds to y = K) we have (3 g = 7^ = beyond one loop, and 
to all orders, in accordance with Ref. [|15|. The results of Eqs. (2.5) 



7£ = 7x 

may readily be specialised to the case of supersymmetric QED simply by setting y 



0. 



3. General Non-Abelian Theory 

We now consider a non-abelian theory with gauge group Q and superpotential 




where Xh $ are multiplets transforming under the S*, S and adjoint representations of 
Q respectively. For notational simplicity we take the representation S to be irreducible for 
the time being; we shall present the results for a reducible representation in due course. 
In addition to diagrams similar in form to those computed earlier in the abelian case and 
shown in Fig. 1, the two-point function for the vector superfield includes the additional 
diagrams depicted in Fig. 5, because the <p field now has gauge interactions. There are 
also further diagrams involving the gauge coupling g only, which we shall be able to avoid 
computing. 

The diagrams of Figs. l(f),(g) and Fig. 5 contain no internal vectors, and produce 
contributions to f3 g which contain the Yukawa coupling y (apart from the zero-bubble con- 
tribution to Fig. 5(a)). They are in fact the only graphs which contribute y dependent 
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terms up to 0(1/Nf). Graphs with no vector propagators will be the same in the back- 
ground field gauge as in an ordinary gauge. Now in the background field gauge, (3 g is given 
by the vector-field two-point function even in the non-abelian case. We deduce that the 
contribution to f3 g at 0(1/Nf) which contains y is correctly given by Figs. 1(f), (g) and 
Fig. 5. Upon performing the .D-algebra, these diagrams all give rise to Feynman integrals 
with bubble sums as depicted in Fig. 3. The contributions from the diagrams of the form 
Figs. l(f),(g) and Fig. 5 are listed in Table 3. 




(a) (b) (c) 

Fig. 5: Additional vector two-point diagrams for the non-abelian theory. 



Diagram 


Bubble Sum 


Group factor 


1(f) 


m 


2gKyC s 


1(g) 


C{y) 


-gKy[C s - \C{G)] 


5(a) 


B{y) 


\gKC{G) 


5(b) 


-B(y) - 2D(y) 


\gKfC{G) 


5(c) 








Table 3: 0(1/Nf) contributions to (3 g in the non-abelian theory. 



The second column of Table 3 requires some explanation. Firstly, we have defined 
y = yT(S), where T(S)5 a b = tr(S a Sb). Secondly, care needs to be taken with the graphs of 
Fig. 5, which contain two bubble-chains. This leads to extra combinatorial factors. In fact, 
upon performing the D-algebra, Fig. 5(a) yields a Feynman integral with the topology of 
Fig. 3(B), but the calculation differs from that in Appendix A, in that Eq. (A. 8) is replaced 
by 



B(k) -> B(k) = ^(n + l)G n K n = ^-[kG(k)]. 



n=0 



(3.2) 



Similarly, Fig. 5(b) yields two Feynman integrals. One again has the topology of Fig. 3(B), 
but now with 



B(k) - B(k) = ;>> + l)G n+2 K n = — — [G(k) - 1] . 

n=0 ^ > 



(3.3) 
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The other is depicted in Fig. 3(D), and yields a bubble-sum contribution given in 
Eq. (A.25). In the third column, we have defined C(S) = C5I, where C(S) = S a S a , 
and used tr ^g) - = Cs for an irreducible representation (where r is the number of group 
generators). Adding the contributions in Table 3, we find 



Hf Jo M f 



G(x) dx 



C(G). 



(3.4) 



As argued earlier, this represents the full y-dependent contribution at 0(1/Nf), together 
with a single y-independent term which represents the (one-loop) contribution from the 
zero-bubble part of Fig. 5(a). We can then infer the full non-abelian result by using the 
afore-mentioned fact that there are no divergences beyond one loop for J\f = 2, i.e. when 
we set y = K, together with the fact that the part of (3 g we have not yet calculated is 
purely a function of g. We also need to include the rest of the one-loop j3 g . The result is 



P g =gK 



+ 



gK 
N f 



T(S) + ^ / (1 - 2x)G{x) dx 

Nf Jy 

/ G(x)dx-1 C(G), 

J V 



(3.5) 



where K = KT(S). 

The chiral superfield anomalous dimensions are given by diagrams of the form of Fig. 4 
together with the additional diagrams of Fig. 6. 





(b) 



(c) 



-CD- 

(d) 

Fig. 6: The additional Feynman diagrams for the anomalous dimension 
in the non-abelian theory. 
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The individual contributions at 0(1/Nf) are given in Table 4. 



Diagram 


Bubble Sum 


Group factor 


4(a) 


B{y) 


yC(S) 


4(b) 


B(k) 


-KC(S) 


4(c) 







4(e) 


m 


2yyC s 


4(f) 


j,[ B {k)-i]-A{k) 


2yKC s 


4(g) 


i[B(k)-i] 


-2yK[C s - \C{G)] 


6(a) 


B(k) 


-KC{G) 


6(b) 







6(c) 







6(d) 








Table 4-' 0(1/Nf) contributions to anomalous dimensions for the non- 
abelian theory. 

The anomalous dimensions for £ or x are given at 0(1 /Nf) by adding the contributions 
from Figs. 4(a)-(c), while the anomalous dimension for (ft is given up to 0(l/iV/) by adding 
the 0(1/Nf) contributions from Figs. 4(e)-(g), and Figs. 6(a)-(c), together with the 0(1) 
contribution from Fig. 4(d). We obtain: 

1 



7? =7x 



= y + 



N f 

2yC s 
N f 



yG{y)-KG(K) C(S) 



(3.6a) 



G(K) - G(y) + 2 / G(x) dx 



K 



+ ±( y -K)G(k)C(G)-^-C(G). 
N f N f 



(3.66) 



The above results contain as special cases all those presented in previous sections. Once 



again one can check compatibility with the three and four- loop calculations from Ref. [12 



and Ref. f 13 



For completeness we present results for the case when x-i £ transform according to a 



reducible representation. Suppose there are n a pairs of multiplets £| 



each trans- 



forming according to a representation S a , and with ^ a n a = Nf. The superpotential 
becomes: 

xata. T ( 



w 



EE 



ita l got c 



(3.7) 
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where a labels the irreducible representations. It is convenient to generalise the definitions 
of y and K as follows: 



and introduce 



Our results are: 



K 



9' 



8n 2 N 



^2n a T(S a ) 



K 



j a 



(3.8) 



(3.9) 
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K + 2KA^ [ K ( l _ 2x)G ( x)dx _^SJh L (\i-2x)G{x)dx 



N f K Jo 



NfV Jo 



= lXc 



+ w f C(G > 

1 

Nf 



f G(x)dx-1 



y a G(y)-KG(K)\C(S a ), 



i<t> = y + 



2A, 



2A y K 



G(K) -1 + 2 G(x) dx 



K 



(3.10) 



+ 



N f y 
Ky 



G(y) -1 + 2 [ V G(x)dx 
Jo 



■C(G) G(K)-1 G(K)C(G). 

N f K L J N f 



K 



4. Supersymmetric QCD 

We can now deduce f3 g and the chiral superfield anomalous dimension for large- Nf 
supersymmetric QCD from the results of the previous section. We consider supersymmetric 
QCD with Nf pairs of chiral superfields & , Xi> transforming under the S*, S representations 
of Q respectively. In principle, to obtain the results for supersymmetric QCD from those 
for our general non-abelian theory, we need to remove the contributions of all diagrams 
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involving the adjoint field (p. For the anomalous dimension of the chiral field, we simply 
need to remove the diagram Fig. 4(a). The result is then 



7£ =7x 



1 



KG{K)C{S). 



(4.1) 



We note that we may obtain this result simply by setting y = in Eq. ( |3.6o| ). In fact, 
the same principle may also be applied in the case of /3 g . One can convince oneself that 
the vector-field two-point diagrams with a ^-propagator which contribute at 0(1/Nf) are 
precisely those of Fig. l(f),(g) and Fig. 5. These all contain y except for the zero-bubble 
contribution to Fig. 5(a). Hence, we can delete the contributions of diagrams with adjoint 
fields simply by setting y = in Eq. ( |3.5| ), provided we also subtract the contribution of 
the zero-bubble contribution to Fig. 5(a) (i.e. the y-independent term in Eq. ( |3.4j) ) by 
hand. The result is 



P 9 =gK 



T(S) + — [ (l-2x)G(x)dx 
N f Jo 



+ 



gK 

N f 



K 



G(x) dx 



(4.2) 



C(G). 



For the case of SU(N C ) with Nf flavours we find 



TV f K 

P 9 =b K ~ WW9 R + 2 9 K W- / ( X - X ^ G ^ dx 
ZN f N f J 



gK 

N c N f Jo 



K 



(4.3) 



(l-2x)G(x)dx. 



Note that in this case K = K/2. 

It is also of interest to consider the case of supersymmetric QCD coupled to a singlet 
4> by a superpotential 

W=^J ) J2&X l - (4-4) 



N f 

The additional diagrams contributing to the vector two-point function are of the form of 
Figs. l(f),(g). Now these diagrams contain no internal vectors, and hence would be the 
same in the background field gauge as in a conventional gauge. But in the background 
gauge, j3 g is determined by the vector two-point diagrams. Hence the additional contribu- 
tion to (3 g can be computed from these diagrams. The contributions can be read off from 
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the corresponding entries in Table 1 by replacing K by K in the last column. Hence {3 g is 
now given by 

2gKC s rK 



P g =gK 



2 f y 

1 - — / (1 -2x)G{x)dx 



N 



+ 



gK 
N 



f 



f Jo 

k 



+ 



N 



f Jo 



[1 — 2x)G(x) dx 



f 3 
/ G(x) dx — — 
Jo ^ 



(4.5) 



C{G). 



The anomalous dimensions are now given by diagrams of the form Fig. 4. The contributions 
from Figs. 4(a)-(c) can be read off from the corresponding entries in Table 2 by replacing 
K by K in the second column and by KC(S) in the third column, giving 

1 



7? =7x 



N f 



yG(y) - KG(K)C(S) 



(4.6) 



The contributions from Figs. 4(e)-(g) can be read off from Table 2 by replacing K by K 
everywhere, giving 



, 2y 



G(K) - G(y) + 2 / G(x) dx 



K 



(4.7) 



For the case of SU(N C ) with Nf flavours we find 



K 



(3 g =\gK - ^gK + 2gK^ I (1 - x)G(x) dx 



gK 

N c N f 



[ (1 - 2x)G{x) dx - ^ f V (l-2x)G{x)dx. 
Jo N f Jo 



5. Infra-red fixed points 

In this section we compare our result for supersymmetric QCD with perturbation 
theory in the gauge coupling. We begin by giving the known results at one through three 
loops: i 



16tt 2 ^ 



i) 



:i6- 2 ) 2 ^ 



'16tt 2 ) 3 /3( 3) 



AN r - 



N r 



Nr. 



-4N n 



Nf b Nj 



T 2 I 9 



1 



+ 



21N 9 

c 



21- c 1 - 7 



N] Nf 



(5.1) 



2 Recall that our gauge coupling is g/WNf 
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For fig 4 ^ we have only a partial result [ 13 



(16tt 2 ) 4 4 4) 



1 



3 atjv 



+ 



,62 x - 100 „ 6« 

■(— + 2k + 8o)N* + — +4a + 



20 



JJV* 



+ 



36(1 + a)N t 



(34 + 12a)iV c - - ±1 i_ - ( 6 + 36a) * ' 



AT? J AT 3 



(5.2) 



where a is an as yet undetermined parameter, and where k = 6^(3). 



It is easy to show[|T6| that (3 g 2 ^ > when (5^ 



less than but near zero, 



0. For (3 g 

it follows that there exists an infra-red fixed point in the evolution of g under renormali- 
sation i. According to Seiberg[|T3 this fixed point in fact exists for the conformal window 
|AT C < Nf < 3N C . Our expansion clearly demands Nf >> N c and so places the theory 
firmly in the weakly coupled infra-red free regime; nevertheless it is tempting to compare 
our result Eq. ( [4.3| ) with perturbation theory for a value of Nf in the conformal window. 



0.5 



-0.5 



S 

S -1.5 

a) 



-2.5 



-3.5 



1 2 3 4 5 6 7 

9 



10 



Fig. 7: Comparison between Eq. (U-3j) (solid line) and the two, three 
and four loop approximations for (3 g (dashed, dot-dashed and dotted lines 
respectively) , for N c = 3 and Nf = 6. 

In Fig. 7 we plot (3 g against g using 2 ■ • • 4-loop perturbation theory and our result 
Eq. ( |4.3[ ) for N c = 3 and Nf = 6. In the case of f3 g ^ we have set a = 0; based on Ref. JTB | 



This holds also for non-supersymmetric QCD [17 
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we would anticipate that \a\ is 0(1), and the evolution is insensitive to the value of a in 
this region. The IR fixed point occurs at g m 8; of course this result is scheme dependent, 
but it is interesting that it is substantially smaller than that obtained in either the two or 
the three loop approximation; the four loop approximation does not give a fixed point at 
all (unless a is large and positive). 



6. Discussion 

It is quite remarkable that the 0(1/Nf) corrections to the SQCD (3 g depend only on 
simple integrals involving G(x). G has a simple pole at x = 3/2 and consequently (3 g has a 
logarithmic singularity at g 2 / (167r 2 ) = 3/2 and a finite radius of convergence in g. We have 
found that for values of Nf, N c corresponding to the conformal window 3N c /2 < Nf < 3iV c , 
there indeed exists an infra-red fixed point in the gauge coupling evolution, though as we 
emphasised this regime is clearly outside the region of strict validity of our approximation. 

If we had some more terms in the 1/Nf expansion, we could try Pade or Borel- 
Pade techniques and continue to small Nf. This might permit us to extend the radius 
of convergence mentioned above, and also investigate more reliably the conformal window 
3N c /2 < Nf < 3N C . Even the (9(l/iV|) contribution presents considerable technical 
problems, however. It is interesting to note that this calculation would suffice to determine 
the unknown parameter a in Eq. (|5.2j ), and consequently complete the derivation of (3^ 



carried out in Ref. [13]. Perhaps the critical methods of Ref M could be extended to 



superfields and facilitate such calculations. It may also be possible to determine a^9\ 



using Pade approximants[2C]. It would be interesting to see whether the finite radius of 
convergence in g which we noted above persists at higher orders in 1/Nf. In fact, it is 
natural to speculate that the 0(1/Nj) term, for instance, would depend on G 2 , or some 
convolution thereof. 

Since our calculations involve contributions from all orders in perturbation theory, 
we should address the question of ambiguities in DREDjn| which can potentially arise at 



higher orders. In our earlier paper on the present topic||, we argued that the particular 
graphs we consider, consisting of bubble chains inserted onto simple lower order graphs, 
are unambiguous; we also speculated that any DRED ambiguities present in the theory as 
a whole should be equivalent to renormalisation scheme ambiguities and as such can be 
subsumed into coupling constant redefinitions. 
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Appendix A. Bubble sums 

In this appendix we give details of our bubble-sum calculations for the diagrams 
of Fig. 3. We do all calculations with zero external momentum, using supersymmetric 
dimensional regularisation (with d = 4 — 2e) and minimal subtraction (DRED)[^0|. By 
performing subtractions at the level of the Feynman integrals we completely separate the 
calculation of the (subtracted) Feynman integrals from the details of the theory under 
consideration. It is convenient to redefine the <i-dimensional integration measure so that 

d d k 1 

K 2 -( P 2 )- e - (A.l) 



k 2 (k — p) 2 e 

The diagrams we will require are shown in Figure 3. Let us consider Fig. 3(B) first of all, 
as this is the simplest integral to evaluate. As explained in the main text, this diagram 
represents a sum over chains of bubbles of arbitrary length. The n-bubble diagram, before 
subtracting subdivergences, produces a contribution 

—G{e) 4^-r[l + (n+ l)e]r[l - (n + l)e]x^ (A.2) 

where 

°W = r(2-e)r ( j:!ffli + e) (A ' 3) 

and x = /U -2 , \i being the regulator mass. The parameter k subsumes any constant factors 
which will recur on a bubble-by-bubble basis, including a factor of (167r 2 ) -1 for each bubble. 
Since this is an (n + l)-loop diagram, there remains an additional factor of (167r 2 ) -1 which 
we have included in the third columns of the Tables in the main text. Similar considerations 
will apply to Figs. (A),(C), and (D). We must now subtract diagrams with counterterm 
insertions corresponding to each divergent subdiagram. For the n-bubble contribution, 
the divergent subdiagrams consist of subsets of n — r + 1 bubbles, where 1 < r < n 
(not necessarily forming a continuous chain). Such a subdiagram yields a counterterm 

of - — n - r +i j after subtracting all its own subdivergences. The remaining (r — 1)- 

bubble diagram gives a contribution as in Eq. ( |A.2| ), but with n replaced by r — 1. Taking 

into account a combinatorial factor of ( n ^ j for the number of (n — r + l)-bubble 

subdiagrams, we find after subtracting all counterterm insertions that the n-bubble (i.e. 
(n + l)-loop) contribution to this diagram is given by the expression: 

Bn = ^G{e) r-\l ~ re)T{l + re)T(l -re)!^) {-iy +1 x r \ (A.4) 

r=l 
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(We have also incorporated into this expression a factor of (— for an L-loop diagram, 
which derives from the D-algebra.) 



We now write 11 



(l-re)r(l + re)r(l-re)x re = ^L J (re) J . (A.5) 

3=0 

Substituting in Eq. ( |A.4j) , and using the identity 

n+l , v 

A i = E ri " 1 ( r 2i)(- 1 ) r = when j = 1 ^--- n 

r=l ^ ' 

= — (n + l) -1 when j = 



(A.6) 



(which is proved in Appendix B) we find that the pole terms in B n are given by the 
expression 

n n 



where we have written G(e) = ^G n e n . The identity Eq. ( |A.6|) removes all the non-local 
(i.e. lnx-dependent) counter-terms. Now we want to sum over n. In a /3-function or 
anomalous dimension calculation, the result will be given by the coefficient of the simple 
pole in e in the quantity Y2( n + l)-^n° l6 5 which is easily seen to give 



OC 



B(k) = 2^ G n K n = G(k). (A.8) 

n=0 

Next we consider Fig. 3(A). A complication as compared with the calculation for 
Fig. 3(B) is that there are now two types of divergent subdiagram; in addition to those 
consisting of subsets of bubbles, there is also a divergent (?7.-f-l)-loop subdiagram formed by 
erasing the upper propagator in Fig. 3(A). In fact this subdiagram is precisely the diagram 
Fig. 3(B), and therefore produces a counterterm B^° le as in Eq. ( |A.7| ). We therefore find 
that the n-bubble (or (n + 2)-loop) contribution to the diagram is given by 

An _ ^ 1} r-l ( n ^ Hjef 1 r(l - 6 + re)r(l + e - re)T(l + re) g 



+2 7~2 \r -2 ) I- e r{r -1) T(l - 2e + re 



+ e n+2 n + 1 2^^ e J r(i-e)" 

i=0 v ' 



(A.9) 
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where 

if(e)=r(2-e)G(e), (A.10) 

with G(e) as in Eq. (|A.3|) . As for Fig. 3(B), we have also incorporated a factor (— 1) L+1 = 
(— l) n + 3 generated by the .D-algebra. We may now extend the summation in the second 
term to infinity, at the expense of introducing one extra simple pole term together with 
finite terms. We obtain, up to finite terms, 

_ K n ^ / n \ H[e f i r(l - 6 + re)r(l + 6 - re)r(l + re) 

71 e " +2 ^2 V -2 / l ~ e r(r-l) r(l-2e + re) 

K n 1 H(e) 2 T(l + e) e K n 1 
H — zrnr — —r- -x — Ci n +i 



e n+2 n + l l _ e r(l - e) e n+1 

n+2 f , i \ 



" +2 ! _! / w + 1 \ g(e) 2 r(l - e + re)r(l + e - re)r(l + re) „ 

n+le n + 2 ^ ^ ' V 7 "- 1 / 1-e T(l-2e + re) 

K n 1 ij"(e) 2 r(l + e) e ^ n 1 
+ e n+2 n + 1 — e T(l-e) X ~ Tn+1 +1 ' 

(A.11) 

Rather remarkably, we now notice that the second term simply supplies the r = 1 term in 
the summation of the first term. We therefore obtain 

An = _1 ^_ Y { _ 1)r -i r -i + H(e) 2 r(l - e + re)r(l + e - re)r(l + re) ^ r£ 

K n 1 

C n+ i + finite terms. 

e n + 1 

In the same spirit as for Fig. 3(B), we write 

r(l - e + re)r(l + e - re)T(l + re) 



n+le n + 2 ^ K ' \r-lj 1-e r(l-2e + re) 

r=l v ' ' 



(A.12) 



oc 



E M ^ re ) J - ( A - 13 ) 



r(l-2e + re) 

Substituting in Eq. ( |A.12|) , and using Eq. ( |A.6|) , we find 

K n H(e) 2 1 K n 1 
A - = ^T2" 7 mv ZT^ 1 + finite terms " ^ A - 14 ) 

Clearly from Eq. (|AT3|) we have M = ^(iS)^ » and so > from Ec l s - Q)> we 
find 

A n =— to (1 - 2e)- T7 rG(e) 

e J (n+l)(n + 2) ^ (A15) 

G n+ i + finite terms. 

e n + 1 
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As in the case of Fig. 3(B), the contribution to a /^-function or anomalous dimension will 
be given by the coefficient of the simple pole in e in the quantity XX n + 2)A n , giving 



oo 1 



n=0 



=-E 

n=0 



—I— o 

n^o n+1 



G n +i + 



n + 1 



-G T 



(A.16) 



which may be rewritten in the compact form 



A{k) 



1 



K 



G(k) - 1 + 2 / G(x)cb 



(A.17) 



The calculation of G(/c) parallels very closely that of A(k), and yields 



C(k) = -2k" 1 



G(k)-1+ / (1 + 2x)G(x) dx 



(A.18) 



The computation of Fig. 3(D) is somewhat different, however, and we shall explain it in 
detail. The n-bubble diagram together with its counterterms (in this case simply bubble- 
chains, as for Fig. 3(B)) gives 

n+3 /_-j\r— 1 



r=3 



n 
r — 3 



G(e)r(l-re)r(l + re)x re . (A.19) 



The factor of (n + 1) represents the (n + 1) ways of distributing n bubbles in Fig. 3(D). 
Now we write 

oo 

r(l -re)r(l + re)x re = ^P^rt) 3 . (A.20) 

3=0 

It is easy to see from Eq. ( |A.6|) that 



n+3 , \ 

r=3 V 7 



(A.21) 



for j > 1. This ensures the cancellation of non-local terms and leaves only the j = term. 
Using the result 



n+3 

E 

r=3 



-l) r 



n 
r — 3 



(n+ l)(7i + 2)(n + 3) 



which is proved in Appendix B, together with Pq = 1, we have 



(n + 2)(n + 3) 



(A.22) 



(A.23) 
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The contribution to a /3-function or anomalous dimension is then given by the simple pole 
in J2i n + 3)£>n, yielding 

oo „ 

which may be rewritten 

2 f K 

£)(«) = __/ x G(x)dx. (A.25) 



o 



Appendix B. Summation identities 



In this appendix we prove the identities Eqs. (|A.6|) and ( [A. 22 ) , which play a crucial role 



in disposing of non-local contributions and deriving the final bubble-sum results. Eq. ( |A.6| ) 
may be proved starting from 

(l-aO n = E(-l) P * r ("V 

> — n \ / 



Differentiating Eq. ([B.1D / times (where Z < n — 1), and setting x = 1, we obtain 



= £(-l)V(r-l)...(r-/ + l)(^Y (B.2) 

r=0 ^ ' 

Clearly, by taking linear combinations of Eq. (|B.2| ) with different values of /, we may obtain 

= E(-1)V- 1 (^ (B.3) 
for 1 < j < n. On the other hand, by integrating Eq. ( |B . 1| ) we obtain 

-±_ = £(-!)> + !)-'(;). (B.4) 



Upon changing the summation range, Eqs. (|B.3|) and (|B.4j) are clearly equivalent to 
Eq. (|A.6| ), namely 

n+1 x x 

A 7 - = r 7-1 f n i (-l) r = when j = l,2,---n 

ti \r-lj y J (B.5) 

= — (n + when j = 0. 
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Upon integrating the identity 



x 2 \ 



l-s)» = X>l) P ("W 2 , (B.6) 



r=0 

and then setting x = 1, we obtain 

^r + 3\rj (n + l)(n + 2)(n + 3)' l ' J 
which is equivalent to Eq. ( |A.22| ). 
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